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SUMMARY 


It is shown that the strea mli nes in an angular nei^horhood of the 
surface of an unyawed circular cone with an attached shock wave are, 
to a first approximation, portions of hyperbolas - This fact is used 
as a basis for the development of an approximate solution in which the 
shock-ira-ve orientation and the flow field behind the shock wave are 
given explicitly in terms of the free-stream Mach number, the vertex 
angle of the body cone, and the ratio of specific heats of the gas. 

The approximate solution is compared with other approximate solutions 
for the cone. 


INTRODUCTION 


It is well known that within calculable limits the adiabatic steady 
flow of an inviscid perfect gas over a nonlifting circular cone assianes 
a relatively singjle form. The essential characteristic of this type of 
flow is an attached cone-shaped shock wave behind which the velocity and 
state of the gas are constant on conical surfaces whose axes and ver- 
tices coincide with those of the body cone and the shock-wave cone. 

The basic features of this type of supersonic flow, together with 
a graphical method for calculating the flow field by starting with a con- 
ical shock wave of given vertex angle and entering air at a given ^fe.ch 
number, were described by Busemann (ref. l) . Taylor and ^fa.ccoll (ref. 2 ) 
conducted a ccm^jrehenslve investigation of the problem in which exact 
equations governing the flow were developed and solved nimierlcally. 

The mmierical results were shown to bear excellent agreement with the 
experimentally determined values. Subsequently, MaccoU (ref. 5) 
extended the numerical conputations to cones with larger vertex angles 
by using solutions in series for the velocity conponents in the flow 
field behind the shock wave and calculated a maximum possible value for 
the vertex angle of a cone with an attached shock wave. The most exten- 
sive tables and graphs of the flonv- variables, which were prepared on the 
basis of the work of Taylor and ^fe.ccoll, are contained in reference 4. 
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Even thoTi^ axially symmetric flow over a cone with an attached 
shock wave is one of the few compressible-flow problems which admit of 
precise solutions with reasonable congjutational labor, approximate 
solutions in closed form sire of use in certain cases. Approximate 
solutions for the cone may be used in the rapid calculation of super- 
sonic flow over pointed bodies of revolution. As a second example, an 
approximate solution relating the velocity field behind the shock wave 
to the velocity at the body-cone surface may, in some instances, be 
combined with information on strong shock waves in a real gas to obtain 
estimates of the effects of vibrational excitation and dissociation on 
the flow over conical tips in the hypersonic range. 

The first approximate solution for the cone resulted from the line- 
arized treatment by Von Karm^Q and Moore (ref. 5) of the axially sym- 
metric supersonic flow over bodies of revolution. In addition to the 
further development of linearized theoiy (for exaniple, ref. 6 by 
T.i gh tbi 11 ) , Van I^ke (ref. 7) lias develop^ a second-order theory in 
which a first-order solution of a given problem constitutes the starting 
point for an iteration procedure aimed at obtaining closer approxima- 
tions. As one mi gbt expect, the increase in accuracy is acccmpanled 
by a considerable Increase in computational effort. 

Since the linear and second-order theories are based upon the 
assnmpt.inn of small . perturbations and flow that is free of shock waves, 
they are definitely limited, particularly in the realm of hl^ Mach 
numbers, both in accviracy and in range of applicability. 

The conlcal-shock-expanslon method of Eggers and Savin (ref. 8) for 
calculating the flow over pointed bodies of revolution at hi^ super- 
sonic airspeeds utilizes approximate solutions for the cone to deter- 
mine the nature of the flow in the nei^borhood of the vertex. Since 
the existence of an attached shock wave is essential to these approxi- 
mations, they are not subject to the same limitations as the perturba- 
tion methods and, in fact, achieve their best accuracy in the hi^ 

Mach number range. Of the three approximate solutions presented for 
the cone, only the third (ref. 8, appendix A) can be applied with ease 
comparable to that of small-perturbation solution of Von 
and Moore. 

In the present report, it is shown that in an angular neighborhood 
of the body cone the streamlines are, to a first approximation, por- 
tions of hj^rbolas. In keeping with the conical nature of the f lew- 
field behind the shock -wave, the hyperbolas of this family are magni- 
fications of one another with respect to their common center, the cone 
vertex. The generators of the body cone (-the whole lines, rather -than 
•the half -lines which terminate at the vertex) are the asyn 5 >to-tes of 
the hyperbolas of -the family. Thus, the conjugate axes of the hyper- 
bolas all lie on the cone axis, and the ciurves themselves are symmetric 
about the plane normal to the cone axis at the -vertex. 
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In the development of the present approximate solution for the non- 
lifting cone with an attached shock wave, the streamlines behind the 
shock wave are assumed to he hyperbolas of the form just described. 
Without further approximation being necessary, it is found that an 
approximate solution can be derived which yields expressions for the 
shock-wave angle, the pressxire ratio across the shock wave, and the 
surface velocity in terms of the free-stream ffech number, the cone 
angle, and the ratio of specific heats. 


SYMBOLS 


a local speed of sound, referred to maximum velocity c 

c maximum or limiting velocity which would be attained by 

expanding adiabatically into a vacuum 

Cp pressure coefficient 

M free-stream Mach number 

p static pressiare 

u radial velocity con^onent (as in polar-coordinate system), 

positive outward along ray from cone vertex, referred to 
maximum velocity c (fig, l) 

V normal velocity component (as in polar-coordinate system) , 

positive with increasing 9, referred to maximum velocity c 
(fig. 1) 

Vj, reidial velocity con^ionent normal to cone axis (as in axially 

symmetric cylindrical coordinate system) , positive outward 
frcM cone axis, referred to TnayTmum velocity c (fig. l) 

velocity component parallel to cone axis, positive in positive 
direction of free stream, referred to maximum velocity c 
(fig. 1) 

x,r axially symmetric cylindrical coordinates with origin at cone 

vertex and x-axis coincident ^ri.th cone axis (fig. l) 

7 ratio of specific heat at constant pressure to specific heat 

at constant voliome 

e^,e^ third-order remainder terms, as described in text 
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0 angle 'wMch ray emanating iircm cone vertex makes with cone 

axis (fig. 1) 

T s tan 9 
Suhscripts: 

s refers to conditions at hody-cone surface 

w refers to conditions immediately behind shock wave 

1 refers to free-stream conditions (not used in case of free- 

stream J&ch number) 

APPROXIMATIOW TO TBE FORM OF THE 
STREAMLINES NEAR THE SURFACE 


The second-order ordinary differential equation which Taylor and 
JfeccoU derived as governing the conical flow field about a nonlifting 
circular cone with an attached shock wa,ve is equivalent to the following 
system of first-order dif ferential equations (see, for example, ref. 4 , 
pp, vii and vlii, or ref. 2, pp. 28l and 282): 



2 2 

dv _ uv - a ( 2 u + V cot 9 ) 
d0 - v^ 


( 1 ) 


where the local speed of soimd a is related to the radial and normal 
velocity conponents u and v (fig. 1) by 


a 


2 




u 



( 2 ) 


Here, for the sake of sinpllcity, the unit of velocity has been taken ' 
as the maximum velocity c which woiild be attained by allowing the 
gas to flcrt'T adlabatically into a vacuum. 

For the purpose of the present analysis, it is convenient to trans- 
form the system of equations (l). Let v^ be the velocity component 
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parallel to tlie axis of the cone (the line whose equation is 0 = 0) 
and positive in the positive direction of the free stream. Moreover, 
let Vp denote the radial velocity component perpendicular to the cone 
axis and positive outward from this axis (fig. l) . Since the flo^f is 
conical, v^ and are functions of 0 . The dependent variables v^ 

and Vj, axe related to the dependent variables u and v by the trans- 
formation 


u = Vjj. cos 0 + Vj, sin 0 


( 3 a) 


V = -v^ sin 0 + Vj, cos 0 


whose inverse is 


v^ = u cos 0 - V sin 0 


( 5 b) 

(hsi) 


Vj. = u sin 0 + V cos 0 


(4b) 


With the change of dependent variables defined by equations (3) , 
the system of differential equations (l) takes the form 


. dVr 

tan 0 ^ - 

and equation (2) becomes 


a^. 


£_ 


- ^Vj. cos 0 - Vjj sin 0)' 


(5) 


a 


2 


7-1 

2 



(6) 


Dividing the system of eqiiations (5) by 1 + tan ^0 and noting that 

(1 + tan^0)d0 = sec^0 d0 = d(tan 0) 
gives, in conical coordinates, 

, T (■ 

where t = tan 0 . (The fact that this form of the equations has also 
been given in ref. 9 ^ P* 556 , was recently noted.) The system of 
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differential equations ( 7 ) is, then, eqtil valent to the second-order 
differential equation derived hy Taylor and ^fe.ccoll. 

Even though the system of equations ( 7 ) apparently does not admit 
of a solution in closed fonn, certq.in interesting and useful facts about 
the existence and nature of a solution are obtainable from the theory 
of functions of real variables. Consideration is restricted to the 
cases of physical interest, that is, 0 < Tg < 00, 1 < 7 < 00, and 

0 < Ug < 1. 

It follows from the theory of functions of real variables that a 
unique solution of eqmtions ( 7 ) exists in a nei^borhood of t = Tg 
such that the functions and Vj. and their first derivatives with 

respect to t are continuous functions. (See, for exanq)le, ref. 10, 

P* 557*) Therefore, by successively differentiating equations (7), it 
is evident that the hi gh er derivatives are also continuous in a nel^- 
borhood of Tg. These results in^jly that the velocity ratio Vj./vx 
and its derivatives with respect to cot 0 = 1/t must be continuous 
in a neighborhood of Tg. Consequently, it is possible to apply Taylor’s 
theorem with a remainder (ref. 10, p. I05) to the development of the 

function v^/vx la powers of cot 9 - cot 0s = ^ ~ evaluation 

of the coefficients for this development is a matter of strai^tforward 
calculation utilizing equations ( 7 ) together with the boundary condition 


vr,s 

''x,s “ 


(8) 


which is equivalent to the physical statement that the flow at the cone 
surface is along a generator. In particular, the term of zero order in 
the Taylor development is given directly by equation (8) . The coeffi- 
cient of the first-order term may be determined as follows: 


d 

d(l/T) 



/_t2 a . Zr \ 

I Wg 



^r 



Vj.2 2 . d.Vy 

„ 2 dT 
X 
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Similarly, the second-order coefficient is found, after seme calculation, 
to he 


1 

d^ Vr 


2 J 

d(l/r )2 vx 

B 1 


The Taylor development, therefore, has the form 



(9) 


where is of the order of the third power of the argumentj that is. 



Eq.iiation (9) can he written as 



1 + 


ts/T 

1 + Tg2 


(t 



+ £5 


(9a) 


where is also of third order. In cases in wtiich the difference 

T - Tg remains small throughout the flow field behind the shock wave 
(that is, where the shock wave lies close to the cone surface), equa- 
tion (9a) can he approximated hy 



( 10 ) 


This relation constitutes the basis of the approximate solution to be 
developed. Of secondary in^iortance in their effects upon the accuracy 
of the approximation expressed in equation (lO) are (a) the tendency, 
most pronoimced in the case of a cone of small vertex angle, of the 
factor Tg/x partially to offset the increase of error with increasing 
and (b) the slight tendency of the factor l/(l + ’’^s^) improve the 
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acciiracy with increasing cone angle. The effect of 7 enters only in 
terms of higher order. (See, for example, ref. 5^ P- ^ 66 .) 

The results of the foregoing analysis may he stated as follows: 

In an angular nei^horhood of the hody-cone surface, the variation of 
the velocity ratio v^J/v^ with the variable t = tan 0 is, to a first 
approximation, given by equation (lO) . 

Consider now an axially symmetric cylindrical coordinate system 
■:-d.th origin at the vertex of the body cone suQh that the x-axis coin- > 
cides with the cone axis and has its positive direction, downstream, and 
let r denote the distance of the point (x,r) in this system from the 
axis. Then, 


r 

— = T 
X 


( 11 ) 


In this coordinate system, the differential equation of a streamline is 

^ ^ Vp 
dx Vx 

With the use of eqmtion (U) and the approximate equation (lO), the 
differential equation of an arbitrary streamline becomes 

dr _ _ 2 X 

^ T 

Separating the variables and integrating gives, as the equation of the 
family of streamlines, 

r^ = Tg^x^ + Constant (12) 

Men the arbitrary constant of integration is permitted to take on 
all positive values, equation ( 12 ) represents a family of geometrically 
similar hyperboi^s with centers at the origin (the cone vertex), trans- 
verse axes throu^ the origin and normal to the cone axis, and common 
eccentricity esc 9s- The hyperbolas of the family have common asymp- 
totes which are generators of the body-cone siarface and whose equations 
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For the set of positive values of the integration constant, this family 
of curves completely fills the portion of space exterior to the solid 
hody cone and its natural counterpart which is obtained hy extending 
the cone generators throu^ the vertex. The conjugate hyperbolas, 
obtained hy assigning negative values to the arbitrary constant in equa- 
tion (12) , lie inside the cone and are irf.thout physical significance in 
this case. The portions of these hyperbolas upstream of the shock wave 
are discarded and replaced hy parallel flow in the subsequent formulation 
of the approximate solution for the flow over cones. 

In the preceeding analysis, it has been shown that in an angular 
nel^horhood of the cone surface the streamlines are, to a first approxi- 
mation, portions of geometrically similar hyperbolas whose common cen- 
ter is the cone vertex and whose asymptotes are extended generators of 
the cone surface. 

The geometric form of the streamlines near the surface depends, to 
a first approximation, solely on the form of the surface (that is, the 
vertex angle of the cone) , but is independent of the velocity at the 
surface and, therefore, independent of the free-stream Ifech number. It 
has been pointed out to the author that equation (12) can also be obtained 
by applying the one -dimensional continuity equation to a stream tube 
adjacent to the surface of the body cone. In this connection it is noted 

that, to a first approximation, the resultant speed ^/u^ + v^ = + v^.^ 

is constant near the body-cone sTjrfacej since the flow is isentrqpic, a 
similar statement applies to the density. However, it should be noted, 
also, that constant resiiltant speed and constant density are not a part 
of the approximate solution to be developed. 


DERIVATICfN OF APPROJOmTE SOLUTION 


In order to develop an approximate solution for axially symmetric 
flow over a cone with an attached shock wave, eqiiation (lO) is now 
assumed to hold throughout the flow field behind the shock wave. If the 
coefficient t in the first of equations (7), namely, 


dr 



is replaced by its value as determined by equation (lO), the resulting 
differential equation can be written in the form 



^x + -^s 


loge Vj.j = 


0 
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Hence, 



= Constant 


Thus, 



( 15 ) 


When the approximate equations (lO) and (15) are solved simultan- 
eously i-rLth the use of equation (8), the resulting approximate equations 
for the velocity components v^ and v^. can be t/ritten as 


Vx,s 


( tan 9 
I tan 0s 


sin^0s 


(l4a) 


tan 0 


, tan 0 


2 

cos 0g 


(I 4 b) 


Since the boimdary values of the velocity components are related to the 
res^lltant velocity Ug at the body-cone surface by the equations 


’^x,s = "s cos 0s 


(15a) 


■^r,s = sill 0g (15b) 

it is seen that equations (l 4 ) essentially give the velocity components 
in terms of the variable 0 and any two of the four parameters s^ 

Vr,s^ Ug, and 0 g. (in the solutions of refs. 2 , 5^ und 4 , the param- 
eters which have been chosen are Us and 0 s-) 

In deriving the approximate relations (l 4 ), use was made of the 
first of equations ( 7 )j but not the second. The significance of neg- 
lecting the second of equations (7) is discussed in appendix A. 

The approximate equations (l 4 ) relating the velocity field behind 
the shock wave to the velocity at the body-cone surface may, in some 
cases, be combined with information on strong shock waves in a real gas 
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to obtain estine,tes of the effects of vibrational excitation and dis- 
sociation on the flow over conical tips in the hypersonic range. Such 
estimates will be most accurate when the relaxation distance measured 
along a streamline is everywhere negligibly small in comparison with 
the length of the conical tip measinred along a generator, that is, when 
the distortion of the conical form of the flow field caused by relaxa- 
tion effects is insignificant. 

In returning to the developnent of the theory based upon, the 
approximate equa-Eibn (lO), consideration is again restricted to perfect 
gases j thus, in particular, the change in state of the gas across the 
shock wave is assumed to be of the Rankine-Hugoniot type and the effects 
of vibrational excitation and dissociation are assumed not to occmr. 

The theory of the oblique shock wave gives, for the velocity conqjonents 
^x,w "'^r,w immediately behind the wave, the exact eqmtions 



In line with the previoiis notation, the velocity components 

Vj, ^ are parallel and peri}endicular, respectively, to the free-stream 

dlrectlonj as before, they are referred to the limiting velocity c. 
The free-stream tfech number and the shock-wave angle are denoted by M 
and 9^, respectively. 

Since equation (lO) is assumed to hold throu^out the flow field 
behind the shock wave, immediately behind the wave it becomes 

~^x,w ^ ^ 1 

Vj.,w tan 9^ ^^2 sln^Q^ 


1 
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T-Jhen the velocity ccm 5 )onents in this equation ane replaced hy their 
values as given in equations (l 6 ), the resulting approximate equation 
is 


2 

7+1 


^sln^e^ 



sin^Gg 


or 


sln^0w = ^ + ^ ^ sinSs (l?) 

if 2 

Tn this approximate equation the shock-wave angle^ the free-stream Mach 
nTjmber, the hody-cone angle, and the ratio of specific heats are related 
in such a way that the equation may he solved for aiDy one of the four 
parameters in tenns of the other three- It is evident that the equation 
becomes exact in the limit as the hody-cone angle approaches zero. 

An interesting theoretical comparison between the cone and the wedge 
can be obtained by Tising the approximate equation (l?) i'ts exact 
cq\interpeLrt from the oblique-shock-wave theory (ref. 11 , p. 57 )* For 
infinite Mach number, it foiLlows from equation (l?) that the limit of 
the ratio of the shock-wave angle to the cone angle, as the cone angle 
approaches zero, is given by 



Thp corresponding limit in the case of the wedge is easily found to 
be (7 + D/2. 

T+iP approximate eqmtion (I 7 ) , for given values of the free-stream 
Mach number M and the cone angle 83 , gives essentially one signifi- 
cant value for the shock-wave angle 0^. Moreover, this value is, in 

general, a reasonably good approximation to the physically observed or 
"first solution" value of 9w, but not to the larger "second solution" 
value. (See ref. it.) This is consistent with the limitation of the 
applicability of equation (lO) to cases in which the shock wave lies 
near the cone surface. 

The approximate expressions that result for the velocity coniponents 
immediately behind the shock wave, when equation (l?) is used to replace 
Sv in equations ( I 6 ) , are 
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cos 0s 


(I8a) 



Eqiiations ( 18 ), then, give the approximate values of the velocity com- 
ponents immediately behind the wave in terms of the cone angle and the 
free -stream lyfe-ch nimiber. 


Using equation (l4a) , evaluated jvist behind the shock wave, together 
with equation (l5a) yields the approximate relation 


Us = 


'x,w 


'tan 


cos 0R \tan 0, 


( 19 ) 


wy 


\Jhen equations ( 17 ) and (l8a) are used to eliminate &w and Vx,w from 
equation (l9)^ the resulting expression for Us is 

Bin^e^ 


1 2 
— M 


1/2 


^s = 


[1 + 2LJL_1 

2 


(cos 0s) 


M - I 1 + 


7 + 1 2 
• xr 


M sin 0s 


1 + 2U_1 i^Bin^0. 


tan^08 


( 20 ) 


In this approximate equation for the velocity at the surface of the 
body cone, as in the case of equations (I 8 ), the parameters are the 
free-stream Mach number M and the cone (semi vertex) angle 0g. 


Thus, equations (l^l-), (15), (IT)^ (I 8 ), and (20) give an approxi- 
mate description of the velocity field and the orientation of the 
attached shock wave in terms of the cone angle, the free-steam Mach 
number, and the adiabatic constant of the gas. 
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The pressure ratio acaross the shock wave is given hy the exact 
equation (ref. 11 , p. 57 ) 


•n_ 7 + 1 W 


7-1 

7+1 


where p^^ is the static pressure in the free stream ahead of the wave 
and p^ is the static pressure immediately behind the wave. Replacing 

sin^G^ in this equation hy its -approximate value as given in equa- 
tion (17) gives the approximate expression 


^ = 1 + 7l^sin^0g (21) 


Since the flow from immediately behind the shock wave to the surface of 
the body is assumed to be isentropic. 



1 - 


1 “ (Tx,w) - (‘'^r,w) 


7 

7-1 


( 22 ) 


where Pg is the static pressure at the surface. 

Consequently, the surface pressure coefficient Cp,s^ which is 
given by 



can be calculated approximately, with the use of equations (18), (20) , 

(21) , and (22) , for given values of the body-cone angle, the free-stream 
Mach number, and the ratio of specific heats. 

Althou^ it is apparently not worthwhile in the general case to 
reduce this set of eq\iations (that is, equations (18) and (20) to (25)) 
to a single expression for the pressure coefficient, considerable sim- 
plification does result when the free-stream Mach number becomes infinite. 
This special case, which is of some theoretical interest, is discussed in 
appendix B. 
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COMPARISON WITH OTHER SOLUTIONS 


In order to compare the present approximte solution with the 
three approximate solutions of Eggers and Savin C^ef . 8 ) , computations 
have been canried out for air with 7 = l-405 in a number of cases 
which represent rather wide ranges of the cone angle and the free-stream 
ffech number; the cases have been selected so as to correspond to exact 
solutions of reference 4. The approximate methods were used to calcu- 
late the shock-wave angle the siirface velocity Ug, and the sur- 

face pressiare coefficient Cp^s* 


The first approximate solution of Eggers and Savin for the cone is 
given in the form of equations (12) and (I 9 ) of reference 8. The first 
of these equations is an approximate expression for the flow-direction 
angle. The ^fe.ch number parameter occvirrlng in the equation is evaluated 
by solving the equation simultaneously with the appropriate oblique- 
shock-wave equations; a semlgraphlcal method, or its equivalent, may be 
used to do this. Wh^ the necessary parameters have been determined, 
the tvro equations give the flow direction and speed, respectively, in 
the flow field behind the shock wave as a function of the position 
angle 0. The solution is primarily applicable to slender cones. 

Eqiiations ( 27 ) and ( 30 ) of reference 8 make up the second approxi- 
mate solution presented by Eggers and Savin. The ccunputational detalln 
are similar to those of the first method. The second solution is pri- 
marily applicable to cases in which the shock wave lies close to the 
body-cone surface and, in this respect, the solution is similar to the 
present approximation. Reference 8 indicates that the second method, 
when it is applicable, is to be preferred to the first method. 

The third approximate method given by Eggers and Savin (appendix A 
of ref. 8) corresponds to the limiting case of the solution of the 
present report for small cone and shock-wave angles. 


In a note on the relation of the hypersonic similarity law to eixlally 
symmetric flow over cones, Lees (ref. 12) has developed an approximate 
solution which is comparable to the third method of Eggers and Savin 
but which is, in general, less uniformly accurate in the calculation of 
both the shock- 1 -ra.ve angle and the surface press\ire coefficient. The 
approximate expression obtained by Lees relating the wave angle to the 
cone angle, the free-stream Mach number, and the specific-heats ratio 
can be put into the following form: 



7 + i 

2 . 
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which permits an easy comparison with equation (A 7 ) of reference 8 , that 
is- 



7+1 

2 



and -with equation (17) of the present paper. 

The conqjarison of the three approximate solutions of reference 8 
■with the approximate solution of the present paper is given in figure 2. 

In order to facilitate the con5)arison of the accuracies of the methods, 
the approximate results are presented in the form of ratios of the cal- 
culated -values to the exact -values as given in reference 4 . 

Figure 2 generally sho^vs that, for 0 s = 5 °^ the present method is 
considerably less acciurate -than -the lengthier first and second methods 
of reference 8 hut is approximately equi-valent in accuracy to the third 
method of reference 8. For 9 s = 20 °, the present method is sli^tly 
less accurate than the first and second methods of reference 8, hut is 
considerably more accm-ate than the third method. For 0 g = 55 °^ the 

present method is slightly more accurate than the second method of 
reference 8, ajid is considerably more accurate than the third method. 

Figure 5 shcnvs -that the present approximation conqjares favorably 
in both accuracy and xmifonni-ty of accuracy with the first- and second- 
order or small-pert\arbation solutions. 

For given -values of the cone angle 0 s and the ratio of specific 
heats 7, equation (20) may be used wi-tli successive approximations to 

determine the free-s-tream Mach ntnriber which corresponds to Us = 

that is, to a Ifech nimiber of unity on -the surface of the body cone. 

These calculations have been carried out with 7 = 1.405 for cones 
ranging in semivertex angle from 0° to 50°, and -the results are presen-ted 
in figure 4 together -wi-fch the corresponding exact data (ref. 4 ). The 
two curves, approximate and exact, give, as a function of the cone angle, 
the minimum free-stream tfech number for which the flow is entirely 
supersonic. 

A detailed examination of -fche accuracy of the present approximate 
solution in calculating the shock-wa-ve angle &w, the surface velocity Us, 
and' the surface pressure coefficient Op g is given in figure 5 for 
7 = 1.405. The conq)utations were performed for several cones with 0 g 
ranging from 5° to 50 ° for Mach numbers corresponding to those for which 
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exact calculations were availatle in reference 4. The lower Mach mam- 
her limits of the curves correspond to a surface tfe-ch numiber of unity 
(consistent with the approximate cturve of fig- 4) . The tendency towEird 
increased accuracy at the hi^ Mach numhers is evident in all the cases 
shown in figure 5 - 


For the limiting case of infinite Nfech numher, which is discussed 
in appendix B, the approximate method has heen applied to the calcula- 
tion of the surface pressirre coefficient taken, for convenience, in 
ratio to the Newtonian pressure coefficient (Cp,s) 2 j " ^ sin^Gg. The 
computations were carried throu^ for three values of the ratio of 
specific heats (y = li, 1.405, l^J with 0g ranging from 0° to 50°. 

The results are compared in figure 6 with the exact results given in 

reference 4 for 7 = 1 — and 7 = 1.405, and with tvro points obtained 

5 

by the numerical integration of equations ( 7 ) for 7 = 1 —. 


CONCLUDING REMARKS 


It has been shown that the streamlines in an singular nei^borhood 
of the surface of an \anyawed circular cone with an attached shock wave 
are, to a first approximation, portions of hyperbolas. This fact has 
been used as a basis for the development of sin approximate solution in 
which the shock-^mve orientation and the flow field behind the shock 
wave sire given explicitly in terms of the free-stream tfech number, the 
vertex angle of the body cone, and the ratio of specific heats of the 
gas. The approximate solution has been conpared with other approximate 
solutions for the cone. 


Langley Aeronauticsil Laboratory, 

Natioral Advisory Committee for Aeronautics, 
Langley Field, Va., July 1, 1955* 
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APPENDIX A 


SOME MATHEMATICAL ASPECTS OF THE APPEOXIMATION 


It has heen shown that equations (l 4 ) constitute an approximate 
solution of the nonlinear system of differential equations (7) • More- 
over, the first of equations (7) is satisfied identically hy the approxi 
mate expressions for the velocity components v^^ and Vj. given in eqiia 

tions (l 4 ) . Direct calculation from equations (l 4 ) and (15) yields the 
equation 



dr 


1 + T, 


which is identically satisfied hy the approximate expressions for v^ 
and vr- Therefore, the approximate expressions for Vx and Vr make 
up an exact solution of the following system of differential equations; 


dvx _ ^ 

d^ dT " a2(l + Tg2 


(Al) 


The factor a^ has heen inserted in the numerator and denominator of 
the ri^t-hand es^ression in order to permit a better conqoarison with 
equations (7) . It is now evident that the approximation embodied in 
equation (lO) > which 3 jed to the approximate solution (l^l-) , is equiva- 
lent to approximating the system of equations (7) by setting the fac- 
tors 1 + t 2 and Vr - tvx equal to constants, namely their values 
at the body-cone surface, and solving the resxalting linear system of 
differential equations exactly. 

It is noted from equation (jb) that the factor yr - tvx can 
be expressed as 

Vj. - rvyr = V sec 9 = v ^1 + 


Thus, setting the factor v^ - tv^ in the nonlinear expression of 

equations (7) equal to its value at the surface of the body cone (that 
is, zero) is equivalent to neglecting the normal velocity component v 
in conrpaxison with the local speed of sound a. If this approximation 
alone is made in the system of equations (l) , which is eqiil valent to 
the system of equations (7), the resulting system can be written in the 
form of a second-order, linear differential equation, that is. 
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cot 9 ^ + 2u = 0 
d.0 


(A2) 


\7hen this equation is ei^ressed in terms of cos 9 as the' independent 
variable, it is easily recognized .as the Legendre eqioation for n = 1. 
Its general solution can he written in closed form and is (ref. I 5 ) 


u = A cos 9 + B f i cos 9 loge ^ ^ - 1 ] 

\ 2 ^ 1 - cos 9 / 

where A and B are arbitrajcy constants to be evaluated in terms of 
the boundary values. This approximation, which gives good results 
when v2 « a2^ was recognized as a good approximation in the hypersonic 
speed range a nimiber of years ago by several persons, independently, at 
the Langley Aeronautical Laboratory. The approximation associated with 
the Legendre equations has been found to lead to results which are more 
complex than those which have their origin in the approximate equation (lO) . 

The only approximation which has been introduced in developing the 
approximate solution of the present report for the axially symmetric 
' flow over a cone with an attached shock wave is that of equation (lO) . 

This basic approximation has been found to be equivalent to specifying 
a flow pattern in which the streamlines are portions of hyperbolas of 
the type described earlier. Inasmuch as the second of differential 
equations ( 7 ) was not used in deriving the approximate solution, the 
problem was not overdetermined by the introduction of equation (lO) . 

The relation of the second of equations ( 7 ) to the approximate sblution 
is clear from the form of the approximate system of differential equa- 
tions (A1) which was found to be equivalent to introducing the approxi- 
mation embodied in equation (lO) . 
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APPENDIX B 

THE LIMITING CASE OF INFINITE MACH NUMBER 


For the limiting case of infinite free-stream MacH n\miber, equa- 
tions (18) and (20) to (25) yield the approximate expression 

7 


(S,sj 


M=co _ 


(^p,s) 


N 


1 - cos'^Bp 


7+1 


, 7 +■ 1 • \8iB^es 

1 - sin 0g' 

7+1 2r> 

cos^9g 


27 


sin^Ss 


7-1 


(Bl) 


where 


(S,s) 


'N 


denotes the surface pressure coefficient according to 


the inelastic-impact theory of Newton, that is. 


(S,s)„= 2 sln^e^ 


^•Then sin20s« 1, the application of equation (Bl) is, under normal 
observances of accuracy, subject to rather large coniputational errors. 
For such cases, it is preferable to use an approximate form of equa- 
tion (Bl) obtained by expanding the quantity 


1 - COS^0f 


7+1 2 \slii 0s 

1 - l-r — sin'=^0s 

7 + 1 p_ 

cos'^Of 


in powers of sin^0g. The resulting series is 


^ 7 + 1\ p 

^1 + loge — ^ — J sin'^08 + 


7 - 1 


7+1 1 / 7 + 1\^ [ a 2 ^ \ 

- loge — 2 ^loge 2 / 


Consequently, for sin^0g « 1, equation (Bl) takes the approximate form 
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(S,s)M=a. 



|l + logg 


2 



- loge 


y + 1 



(B2) 


For 0g = 10°, this equation is generally to he preferred to equation (B1) 
in performing numerical computations. 


Finally, equation (B1) permits one to write 


Urn loge 

7—^1 




after the indeterminate form concerned has been resolved. Thus, it 
follows that 


lim 
7 — 3^1 




(B3) 


Since the inelastic-impact theory of Newton is applicable to the flow over 
an unyawed circular cone of a gas -vrith an infinite number of internal 
degrees of freedom, equation (B3) shows that equation (B1) becomes exact 
in the limit as 7 — ^1. Similarly, equation (17) predicts correctly that, 
for infinite Mach number and 7=1, the shock wave lies on the body-cone 
surface. For given values of the velocity and specific enthalpy of the 
gas ahead of the shock wave, the free-stream Mach number necessarily 
becomes infinite in the limit as the ratio of specific heats approaches 
unity. At the same time straightforward analysis more generally shows 
that the approximate expression (17) for the shock-wave angle and the 
approximate expression for the surface pressure coefficient which can be . 
obtained from equations (18) and (20) to (23) both become exact in the 
limit as 7 — >1; that is, they reduce to the corresponding expressions 
in the inelastic-impact theory of Newton. Therefore, the present approxi- 
mate solution for the cone becomes exact in the limit as 7 — >-1 whether 
the free-stream Mach number is restricted to the value infinity or is 
not restricted. 
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Figure 1.- Axially symmetric flov over cone vith attached shock wave 
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(a) Shock-wave angle. 


Figure 2.- Con^jarlson of present method with approximate 

reference 8 . 7 = 1 .^ 5 - 
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(c) Surface pressTure coefficient 
Figure 2 .- Concluded. 
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Free-stream Mach number, M 


Figure 3 -- Surface pressure coefficients for 

angles. 7 = 1.405. 
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Figure 4 .- Free-stream Mach numher plotted against cone semivertex angle 

for surface Mach number of 1. 
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